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Abstract 

We study the third order pion-nucleon scattering amplitude obtained from heavy 
baryon chiral perturbation theory inside the Mandelstam triangle. We reconstruct 
the pion~nucleon amplitude in the unphysical region by use of dispersion relations 
and determine the pertinent low-energy constants by a fit to this amplitude. A 
detailed comparison with values obtained from phase shift analysis is given. Our 
analysis leads to a pion-nucleon u-term of o"(0) = 40MeV based on the Karlsruhe 
partial wave analysis. We have repeated the same procedure using the latest solution 
of the VPI group and find a much larger value for cr(0). 
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1 Introduction 



A detailed understanding of elastic pion-nucleon scattering in the low energy region allows 
for precise tests of the chiral QCD dynamics. Recently, this process has been investigated 
to third order in heavy baryon chiral perturbation theory by various groups. At that 
order, one has to deal with tree graphs and one loop diagrams. At second and third 
order, there are in total four and five so-called low-energy constants, respectively, which 
have to be determined by comparison to data. The loop contribution only depends on 
the weak pion decay constant F = 92.4 MeV and the axial-vector coupling = 1.26. 
In ref. particular combinations of threshold and subthreshold parameters were found 
which do not depend on the dimension three low-energy constants (LECs). Also, it was 
shown that the ensuing numerical values of these dimension two LECs, called q, can be 
understood in terms of resonance exchange. In ref.p], the threshold parameters based on 
the Karlsruhe-Helsinki analysis together with the pion-nucleon cr-term were used to fix 
the nine LECs. The results for the q were in good agreement with the ones obtained in 
[H] . The amplitude given in |^ was used later to obtain S- and P-wave phase shifts below 
the first resonance 0. The most systematic study was performed in ref.0], where the 
S- and P-wave partial wave amplitudes from three different analyses were used (in the 
physical region and in the range of the lowest existing data) to fit the LECs. However, 
chiral perturbation theory is expected to yield the most reliable predictions for s (the 
center-of-mass energy squared) and t (the squared invariant four-momentum transfer) 
lying inside the Mandelstam triangle (depicted in fig.|l|), for essentially two reasons. First, 
in this region the scattering amplitude is purely real and it is well known that at a given 
order in the chiral expansion, the real part is in general more precisely determined than 
the corresponding imaginary part (since the latter only starts at one loop order). Second, 
in the interior of the Mandelstam triangle the kinematical variables t and (s — u)/4m 
take their smallest values.^ As this region is unphysical, there is no direct access by 
experimental data. By the use of dispersion relations this problem can be circumvented. 
This is done here. First, using data from phase shift analysis we construct the pion- 
nucleon amplitude inside the Mandelstam triangle. We then use the chiral third order 
amplitude constructed in |^ to determine the 8 LECs under consideration by a best fit as 
described below and compare their numerical values with the ones obtained previously. 
The ninth LEC is fixed by the value of the pion-nucleon coupling constant used in our 
analysis. We expect a more precise determination of the dimension two LECs since to 
the order we are working, we are not sensitive to the 1/m corrections to the dimension 
three LECs (which are known to be important in the determination of the Cj, compare 
e.g. the values obtained in with the ones of refs.[0, ^ Q.). We remark that Hohler |Q 
has stressed that one should also compare directly the analytical form of the pion-nucleon 
amplitude obtained by different means. We refer to his work for a detailed comparison 
between the dispersive and the chiral representation. Also, it is known that in some small 
regions the heavy baryon amplitude converges slowly. This can be traced back to the 



Throughout, we denote by to = 938.27 MeV and AI — 139.57 MeV the nuclcon and the pion mass, 
respectively, and the Mandelstam variables are subject to the constraint s + 1 + u = 2m^ + 2AP. Fur- 
thermore, we work in the isospin limit to^ = to^ and neglect all virtual photon effects. The exception to 
this are the isospin violation effects contributing to the particle masses, i.e. kinematical effects. 
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fact that the strict heavy baryon hmit tends to modify the analytical structure of the 
ttN amplitude. These effects can be dealt with by subtracting from the amplitudes the 
full Born terms, since the latter generate the singularities (such an effect also appears 
in the discussion of the nucleon electromagnetic form factors, see e.g. 05 H; 9)- K is 
important to stress that this subtraction procedure is not arbitrary since the subthreshold 
expansion of the vrA^ invariant amplitudes is usually formulated by subtracting the Born 
terms to avoid their rapid variations in the appropriate kinematical variables. Another 
way of circumventing this problem is to stick to a relativistic formulation of the matter 
fields, as recently proposed in ref. ||T0[ since in that way all strictures from analyticity are 



automatically fulfilled. Here, however, we are concerned with the comparison of the vrA^ 
amplitude obtained in heavy baryon chiral perturbation theory (HBCHPT) since in this 
framework the most precise predictions for the threshold parameters have been obtained 
and we wish to explore the consistency of these calculations at the order they have been 
performed. In contrast to previous investigations, we confine ourselves to the inside of 
the Mandelstam triangle for the reasons mentioned above. 

The manuscript is organized as follows. Some formalism pertaining to elastic pion-nucleon 
scattering pertinent to our investigation is given in sec. |^. In sec. |^ we construct the 
invariant amplitudes inside the Mandelstam triangle by use of dispersion relations. The 
chiral perturbation theory amplitudes are fitted to these dispersive amplitudes in sec. |[ 
Further results on subthreshold parameters and the cj-term are discussed in sec. ^ The 
summary and conclusions follow in sec. 0. 
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2 Formal aspects of elastic pion— nucleon scattering 

Consider elastic pion-nucleon scattering, 

7r"(gi) + iV(pi)- A92) + iV(p2) , (1) 

with 'a, 6' (cartesian) pion isospin indices and gj, Pi the four-momenta of the pions and 
the nucleons, respectively. The scattering amplitude is usually decomposed in terms of 
the four invariant functions A='=(s,t) and i?^(s,t) (where the superscript '±' refers to the 
isoscalar/isovector part), 



T!^%{s,t) = T+v(.,t)5^'^ + T-^(.,t)-[r^r1 , 



A^is,t) + -iq2 + qiri^B^is,t) 



(2) 



with ^/s the cms energy and t = {qi — 52)^ the invariant momentum transfer squared. 
The Aj {i = 1, 2) denote the helicities of the incoming/outgoing nucleon. In what follows, 
we also need the linear combinations related to the physical channels vr+p ir^p and 
■K~p — > 7r~p. These are related to the isospin amplitudes by 



X±(s,t) = ^(x„(s,t)±X+(s,t)) , X = {A,B} 
The invariant amplitudes fulfill the crossing relations 



(3) 



y4+(s,t,M) = y4+(u,t, S) 

B^{s, t, u) 
A~{s, t, u) 



B iu,t,s) , 
—A^{u, t, s) 
-B^{u,t,s) 



(4) 



For heavy baryon chiral perturbation theory, it is more natural to work directly with the 
non-spin-flip and spin-flip amplitudes (7^(0;, t) and h^iu, t), respectively. For a discussion 
of this point, see 0]. Then, the tcN scattering amplitude takes the form 



E + m 
2m 



ba 



g^{uj, t) +ia ■ {q2 x ) h^{uj, t) 
g'^ (cu, t) + ia ■ (g2 x gl ) /i" (t^, t) 



(5) 



with uo the pion cms energy, 



4^ 



(6) 



E = (g*^ + m^)^/^ the nucleon energy and the pion momenta 

(s - - w?f - Am'^M^ 



^2 ^2 _ ^2 _ 2 

(11 = ^2 = Q =(1 



As 



(7) 



The relations between the invariant amphtudes A^, B"^ and g"^, is 



9 {^,'t) = — — (8) 



h^i^, t) = ' rr rr ' 
0104 — 0203 



with the coefficient functions Cj (i = 1, . . . , 4) given by 



— — m C(j + m Vm^ — + lu"^ 



2 (m2 - M2) 

^2 = , J' -lim^M^ -AM^ -m'^t + M'^t-i 



^3 



4 (m2 - M2) 

+mtu; + 4M^cj^-m (t + 4mu;-4 M^) Vm^ - + u;2| 
1 



(m + - M2 + u;2^ (4 m (m + ^^2 _ ^2 + u;2) - 

C4 = T ^ — . (10) 

4 f m2 - M2 + m + cj2 + (m + cu) Vm2 - M2 + 6^2 j 

This form differs from the one in |^ in that it is exphcitly given in terms of 00 and t 
and also it is more directly applicable for considering processes in regions other than the 
physical s-channel domain. Note that these coefficient functions should not be expanded 
in 1/m when one performs the chiral expansion. Finally, we will also use the crossing- 
symmetric variable u = {s — u)/4m, the scattering angle in the center-of-mass system 
in the s-channel, cos 6' = 1 + t/2g^, and sometimes denote the cms energy by W, i.e. 
s = W^. We also need the subthreshold expansion of the invariant amplitudes with the 
pseudovector Born terms subtracted (as indicated by the "bar" ) |Tl| 



X = Ex„„Z/2m+fc^n ^ X = {A±,5±}, (11) 



m,n 



with A; = 1 (0) if the function considered is odd (even) in u. The Taylor-coefficients Xmn 
are the so-called subthreshold parameters. For a more detailed account of the pertinent 
kinematics, we refer to the monograph ||lT . 



3 Construction of the invariant amplitudes inside the 
Mandelstam triangle using dispersion relations 

In the next step, we use dispersion relations to construct the vrA^ amplitude inside the 
Mandelstam triangle. The Mandelstam triangle is the interior of the region bounded by 
the three lines s = [M + m)^, u = (M + m)^ and t = 4M^ in the Mandelstam plane, 
see fig. |l|. Dispersion relations for vrA^ scattering have been studied in great detail in 
axiomatic field theory in the last thirty years [0, |TT[. Using analyticity, unitarity, and 
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crossing symmetry one obtains dispersive relations for the scattering amplitudes of the 
form 

ReA^{u,t) = -i Im A^{u',t)\ ± } du' , (12) 

ReB^{u,t) = ^rimB^iu',t)(-^T-^]du' + f,{ — T 



(13) 



with z/thr = M + t/{4m) and vb = {t — 2M^)/(4m). Note that the nucleon Born terms 
contribute only to B'^iy^t). Here 

is the pion-nucleon couphng constant squared. There is some debate about its actual 
value, some analysis favoring a somewhat smaller value, p ^ 0.075, as discussed below. 
However, to be consistent with refs.|l|, ^, ^, which mostly use the Karlsruhe phase shifts 
and therefore the value given in eq.(|l4[), we also use the large value for Note that in 



addition we consider an analysis favouring a smaller value for p later on. We mention 
that it seems clear that the imaginary part of the amplitudes A~,B^ fall off sufficiently 
fast so that no subtraction is needed. This might also be true for A'^, see e.g. [1^ . 
However, to precisely determine the subthreshold parameters related to A'^, we perform 



a subtraction for the amplitude A~^ |TT 



Re A+(z/,t) = Re A+(i/ = 0,t) + f ——-^^du', (15) 



where the first term on the right-hand-side of eq.(^) is the subtraction function. This 
will be discussed in more detail below. Note also that the infiuence of any high energy 
contribution is damped further in the once subtracted dispersion relation for A~^{i>,t). 

The absorptive parts of the amplitudes in the dispersive integrals are broken up into two 
parts: a) the low energy part with laboratory momenta fciab < 6 GeV and b) the high 
energy part. The first is constructed by partial waves of the Karlsruhe group (KA84) 



We note that there has been considerable criticism of these partial wave amplitudes re- 
cently (see e.g. the proceedings of MENU 97 [jl5|), however, at present no other analysis 
exist which consistently includes data from threshold to the highest available energies (for 
example, the VPI group uses the Karlsruhe analysis for fciab above 2.1 GeV 0). In our nu- 
merical analysis we work with the S- to K-wave approximation of the physical amplitudes 
A±{i',t) and B±{h',t), i.e. (for clarity, we only display the S- and P-wave contributions 
here. The ellipsis stands for the contributions from the higher partial waves.) 

A+{u,t) W + mf 3/2. ^ , „ a,3/2f J W-m 



A-K E + m 



(16) 

A_{iy,t) W + m [,3/2/ X , N Q /3[/3/2/ X m' 

-|/o+ {ly) + 2/0! (z/) - 3cos^[/i4 (z/) - 2/1I (z/)] 



Att 3{E + m) 
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W — m 



_ ^) ^/i'^'(^) + 2/|^'(z.) - fti'iu) - 2fli'{u)] + . . . , (17) 

^ = ^(^{/^^M + 2i;;/^(.) + 3cos^[A^^(.)+2/,^f(.)] 

+ 3(^!_^) {/y(^) + 2/|^'(^) - (^) - 2A?(^)} + . . . , (19) 

yielding the desired amplitudes A'^,B'^ by eq.(|^). To construct e.g. the amplitude 
i?+(z/, t), we need the partial waves with the appropriate quantum numbers from S31 
up to K315 (in the usual notation l2i,2j)- The pertinent partial wave amplitudes fi±{v) 
for total isospin / = 1/2,3/2, pion-nucleon angular momentum I and total angular mo- 
mentum j = / ± 1/2 are given in terms of the phase shifts and inelasticities tiI^{v) 
via 

fL{^) = ^H±(^) exp5L(^) - 1} ■ (20) 




Figure 2: Argand plot of A^^v.t = 0) for fciab < 6 GeV. All units are appropriate powers 
of the pion mass. 



The corresponding Argand diagrams for A'^ and ujB~ in the low energy region are shown 
in figs. 2 and 3, respectively. These reproduce the ones given by Hohler and collaborators 
with sufficient accuracy for our later purposes, compare with the corresponding figures in 
ref. [|16[. For the high energy part of the amplitudes A~ and B~ we assume a reggeized 
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uj Im B-{iy,t = 0) 




Figure 3: Argand plot oi u B {u^t = 0) for kuh < 6 GeV. All units are appropriate powers 
of the pion mass. 



p-meson exchange in the t-channel of the charge exchange process tt p ix^n to be the 
sole component of the amplitude |]1T[, e.g. 



\mB],{u,t) = D^a{t){a{t) + l}[ — ] , (21) 

with Dq = 38.9 mb, c^o = 1 GeV and the Regge slope is parameterized via a(t) = 0.56 + 
1.08 -t 0. 



In the present work we will make use of the amplitudes along the lines i/ = and t = 
2M^/3 (and also t = 0) to map out the interior of the Mandelstam triangle (cf. fig. 1). The 
choice of this particular value of t is justified in the next section. Consequently, B~{0,t) 
and the subtraction function A~^{0,t) will be of special interest to us. We remark that 
A~{0,t) and B~^{0,t) vanish because of crossing symmetry. A~^{0,t) is reconstructed by 
approximating the real and imaginary part of A'^ in the physical region of the s-channel 
by the available partial waves. Inverting eq.(|T^) then yields 74"'"(0, The result is 



shown in fig. ^. We remark that this amplitude fulfills the Adler consistency condition 

A+(u = 0,t = M\ql = 0,q^ = M^) c^— (22) 

m 

within one percent. Since physical pions do not have vanishing four-momentum, one 
expects a deviation from this relation of the order of M^/(47rF)^ ~ 1.5%. B~{0,t) is 



'^^Note that inside the Mandelstam triangle the amplitude is real. 
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obtained from eq.(|T3D. The i/-dependence of the amphtudes is evaluated in a similar 
fashion by keeping t fixed at the value given above. 



A+{u = 0,t) 




-2'0 -15 -10 -5 



Figure 4: The subtraction function A'^^u = 0,t). The solid (dotted) line refers to the 
KA84 (VPl SP99) partial-wave analysis. 



So far, we have used the KA84 phase shift analysis. We are aware that some of the data 
entering that analysis are now considered inconsistent. Consequently, we will also use the 
latest version of the VPl phase shift analysis, called SP99 |T^, for our calculations. We 
note that the VPl group uses dispersion relations only for some of the amplitudes and in 
a certain energy region, different from the Karlsruhe method. Based on the VPl SP99 
phases, we have constructed the amplitudes A^, as described before. We refrain from 
showing the corresponding figures, but we remark that the Adler consistency condition 
is then violated by 10 percent. The corresponding subtraction function is also shown in 
fig.|| by the dotted curve. Of course, one has to account for the smaller value of the pion- 
nucleon coupling constant used by the VPl group, p ~ 0.076. The Adler consistency 
relation is of importance for the precise determination of the amplitudes when the Born 
terms are subtracted. We believe that for making accurate statements on small quantities 
like the pion-nucleon a term, it is mandatory to fulfill such constraints as e.g. the Adler 
condition to the expected accuracy as explained before. 



4 Chiral perturbation theory amplitudes inside the 
Mandelstam triangle 

We now turn to the comparison with the chiral amplitude. For that, it is important to 
study the dependence of the various invariant amplitudes on the pertinent LECs. We give 

'^^We have been informed by Marcello Pavan that the VPl group is presently improving their represen- 
tation of the subtraction function. 
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here the exphcit expressions for the counterterm amphtudes taken from ref.[0] (in terms 









C2UJ r 


9tt = 






C4 t 


9ct = 


2mF2 



-4ciM^ + 2c2CU^ + C3(2M^-t) , (23) 



+ -^(2M2 - t) {di + rfa) + 4 + — ^ 4 

+ 14 — (2^'-2M2 + t)Ji8 (24) 

, , 2iu ^- ^ , 9a - ^ ^ 

^ct = ^(ai4 -dis) + -^^^rfis , (25) 

^ct = ;;^ + ^' (26) 

to facihtate the later discussion. The //-dependence can be worked out using the fact that 
u and z/ are related via 

u = —(2s + t- 2M^ - 2n?) , (27) 
4m 

and using eq.(|]). Note that most of the LECs appear only in one particular invariant 
amplitude, e.g. ci, 02,03 are fixed entirely from varying g~^{uj,t) inside the Mandelstam 
triangle. Only C4 and dig appear in two different amplitudes. In particular, we expect 
that most of the dimension three LECs cannot be pinned down accurately, since they 
only appear in with small prefactors. This expectation is borne out by the results 
given below. Furthermore, to this order in the chiral expansion we are not sensitive to 
the 1/m corrections to the d^. Such terms are known to play an important role in the 
precise determination of the dimension two LECs, compare e.g. the values obtained in 
ref. with the ones in refs.[0, 0, ^. It should be noted that the amplitudes given in |Q 
are valid for u > M,t < 0. Therefore, when working inside the Mandelstam triangle, one 
has to ensure the correct analytic continuation to the unphysical region u < M of the 
complex-valued one-loop contributions of and by use of |l| 

Vl — = —i \l — \ , arcsin x = — + i In (^x + — l) • (28) 

A similar statement holds for the continuation to t > 0. We use the tree and one-loop 
amplitudes explicitly given in @] but refrain from spelling them out here. 

The chiral amplitudes to order O(g^) do not have the correct analytic structure inside 
the Mandelstam triangle. This is in part due to the Taylor expansion of the lowest order 
amplitude in powers of m~^. To avoid this problem, it is necessary to neglect the tree 
contributions and the counter terms ~ gA- This holds in particular for the tree graphs 
with insertions ~ d\g. The latter allow one to use the physical value of the pion-nucleon 
coupling constant (instead of its value in the chiral limit) in the tree graphs, which are 
responsible for the incorrect analytic behaviour. This procedure is equivalent to working 
with the quantities = — g^ /m and so on instead of A^, B"^ in eqs.(|]),(^. This way, 
all the contributions due to rfig drop out and the latter has to be calculated in a different 
way, e.g. by the so-called Goldberger-Treiman discrepancy (i.e the deviation from the 
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Goldberger-Treiman relation) : 

Fn \ 9a ) 

In the present work we compare the chiral amphtudes with their dispersive counterparts 
in a small region inside the Mandelstam triangle. More precisely, we concentrate on two 
regions. The first of these is concentrated around the point i/ = and t = 2M^/3, the 
so-called "ideal point". In the chiral limit, where the pions are massless, the ideal point 
is given by z/ = t = and the constraint s + 1 + u = 2mo has to be fulfilled, with mo the 
nucleon mass in the chiral limit. Since we are considering physical nucleons and physical 
pions, the ideal point is shifted to t ^ 2M^/3 because now s + 1 + u = 2m? + 2M^. We 
consider two fits, which span the following regions within the Mandelstam triangle 



Fit 1 : around = 0, t = 

with {-0.5 <t< 0.5 , 4:5.7M^ <s< 46. TM^} , (30) 
Fit 2 : around z/ = 0, t = 2M^/3 

with {-0.3 <t< 0.7 , 45.35M2 < s < 46.35} , (31) 

This means that for both fits we cover the range of one pion mass squared in the two 
independent Mandelstam variables. A larger range cannot be accounted for by the third 
order chiral amplitudes. We have convinced ourselves that decreasing these ranges by a 
factor of two does not change our results. Note that for the fits along the line u = 0, a 
simplification arises. Due to crossing symmetry and eqs.(H),(y), each of the amplitudes 
g^,h^ depends on only one of the invariant amplitudes A^,B^. E.g. the non-spin-fiip 
amplitude g~^{i' = 0, t) is related to the modified subtraction function y4+(i/ = 0, t) via 

9H0,t) = ^^^'^ A^O,t). (32) 

A fit to each of the dispersive amplitudes g'^{i',t),h^{i',t) at 21 equidistantly separated 
points in the two directions of fixed v and fixed t in the regions given in eq. (|30|) yields the 
following values for the LECs (we give here the averaged values obtained from fit 1 and 
fit 2 together, the individual results are displayed and discussed below): 

ci = (-0.81 ± 0.12) GeV"^ , Ca = (8.43 ± 56.9) GeV"^ , 

C3 = (-4.70 ± 1.16) GeV~^ , C4 = (3.40 ± 0.04) GeV"^ . (33) 

The errors given are the the ones from the two fits added in quadrature. The errors 
for each fit have been determined from the square root of the diagonal elements of the 
error matrix. We note that there is a sizeable variation in the precision with which the 
dimension two LECs are determined. It is expected that the isovector amplitudes are 
more accurately predicted at third order by chiral perturbation theory since the isoscalar 
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amplitudes vanish to leading order in the chiral expansion.^ That is also the reason 
why the LEG C4 is more precisely pinned down than the LECs Ci^2,3- In table |I| we 
compare the values of the dimension two LECs for fit 1 and fit 2 separately with previous 
determinations. In ref. (sub)threshold parameters and the cr-term (i.e. quantities free 



Ci 


Fit 1 


Fit 2 


Ref.| 


1 


Ref. 


i 




Ref.[ 


1 (Fit 1) 


1 

2 
3 
4 


-0.81 ±0.15 
9.35 ±66.7 

-4.69 ± 1.34 
3.40 ± 0.04 


-0.80 ±0.07 
7.52 ±45.0 

-4.70 ±0.95 
3.40 ± 0.04 


-0.93 ±0.10 
3.34 ±0.20 

-5.29 ±0.25 
3.63 ±0.10 


-0.94 ±0.06 
3.20 ±0.10 

-5.40 ±0.06 
3.47 ±0.06 


-1.23 ±0.16 
3.28 ±0.23 

-5.94 ±0.09 
3.47 ±0.05 



Table 1: The dimension two LECs ci^2,3,4 in GeV~^ from our determination based on the 
KA84 phase shifts compared to previous ones using also information from vrA^ scattering 
in the physical region and the cr-term. 



of dimension three LECs) were used, while in ref. [§] the q and dimension three LECs 
were determined from a fit to threshold parameters of the Karlsruhe-Helsinki group . 
Finally, the determination in ref.lQ (fit 1 therein) used only the available low-energy 
KA84 phase shifts for pion laboratory momenta in the range from 40 to 97 MeV. The 
results for ci, C3 and C4 are very consistent for all the different determinations.^ The 
exception is our value for C2, which is essentially an undetermined quantity. There reason 
for that is that in the amplitude , which is small, the LECs Ci and C3 are weighted with 
factors of M^, whereas the term ~ C2 is proportional to uo"^ (to leading order in the 1/m) 
and this quantity is suppressed by a factor of 10 compared to in the region around 
the center of the Mandelstam triangle. As a consequence, we can not make a precise 
prediction for the combination c"*" = C2 ± C3 — 2ci, which can also be deduced from pion 



scattering off deuterium pO| . One could, however, turn the argument around and use the 
value of c"*" determined in ref.[^, c"*" = —0.09 ± 0.37 GeV~^ Neglecting the uncertainties 
in ci and C3, this leads to 

C2 = (2.99 ± 0.37) GeV"^ , (34) 

which is consistent with the other determinations compiled in table |I]. It is also possible 
to pin down the LECs 01^3,4 from the long-range part of the proton-proton interaction 
based on the chiral two-pion exchange potential. This has been done in ref. and leads 
to results consistent with the ones found here and in previous works. The values obtained 



in ref.|2T| are Ci = -0.76(7), C3 = -5.08(28) and C4 = 4.70(70) (all in GeV-^). Finally, 
we remark that the errors we quote (compare eq. ( ^3]) and table |I]) must be handled with 
care: to our knowledge, there is no error analysis available for the KA84 phase shifts 
and elasticities. Thus, it is impossible to give reliable errors for the dispersive amplitudes 



^^Note that this does not neccessarily imply a similar pattern for the accuracy with which the corre- 
sponding counterterm contributions can be pinned down. This argument can also be invalidated in case 
of zero crossings of certain amplitudes. 

#^Note that we come back to the different value for ci obtained here and from scattering data alone M 
in section |^. 
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inside the Mandelstam triangle. In the present work we assume an error of 10 % for the 
amphtudes of interest. In this sense our errors are arbitrary. Stated differently, using 
other errors does not change the central values but the uncertainties deduced from the 
corresponding error matrix. 

We now turn to the dimension three LECs. We have found that if we fit along lines 
of fixed V or fixed t, they can be pinned down with good precision. However, the so 
determined values are not mutually consistent. The reason for this lies in the fact that 
the third order amplitude employed here is not sufficiently precise to describe the small 
contributions from accurately. In table || we compare our results for the di. Only 
the combination (ii4 — di^ comes out consistent for the determinations within the regions 
given by fits 1 and 2, respectively. Note, however, that the overall description of the 
amplitude is fairly poor, so that this determination of d^ — di^ has to be taken with 
some caution (even though the uncertainty in this LEG from the fits is tiny). It appears 
that the 1/m^ corrections with fixed coefficients and the Cj/m terms are more important 
than the genuinely new operators of dimension three.|^ Also shown in table ^ are results 
from previous determinations determined by fitting to threshold parameters and phase 
shifts. It would be valuable to have additional information on these LECs, say from a 
study of elastic pion-nucleon scattering to fourth order in the chiral expansion or from 
the reaction vrA^ ttttN. 



di 


Fit 1 


Fit 2 


Ref.§ 


Ref.§ (Fit 1) 


1+2 

3 
5 

14-15 


3.33 ±0.01 
-152.3 ±0.24 
-0.11 ±0.01 

0.96 ±0.00 


-1.91 ±0.14 
98.16 ±6.90 
-0.20 ±0.01 
0.93 ±0.00 


2.4 ±0.3 
-2.8 ±0.6 

1.4 ±0.3 
-6.1 ±0.6 


3.06 ±0.21 
-3.27 ±0.73 

0.45 ±0.42 
-5.65 ±0.41 



Table 2: The dimension three LECs di in GeV~^ (and combinations thereof) from our 
determination compared to previous ones using information from nN scattering in the 
physical region. 



Using now the SP99 phases as input to construct the amplitudes and redo the fit, we 
find (taking into account the different value for the vrA^ coupling constant used by the 
VPI group) 

[(-1.53 ±0.18) GeV"^] , 
[(3.22 ±0.25) GeV"^] , 
[(-6.20 ± 0.09) GeV~^] , (35) 

[(3.51 ± 0.04) GeV"^] , 
[(-1.53 ±0.18) GeV-2] , 

where we have given the range based on fits 1,2 as described above. The numbers in the 
square brackets are the results of fit 3 of ref.]^] based on the VPI SP98 solution. The 

"^^A similar result for another set of dimension three LECs was recently reported in ref . p2| . 



C2 
C3 
C4 

di4 + C?i5 



(-3.00... -2.96) GeV"^ 
(-9.51... 7.22) GeV"^ 

(-6.63... -6.20) GeV'^ 
(3.39... 3.40) GeV"^ 
(0.92... 0.93) GeV-2 
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uncertainties of the so determined LECs are very similar to the ones obtained based on the 
KA84 partial waves. Again, we note that C2 cannot be determined. Notice in particular 
the very large value for ci, which will be commented on below. The value for C4 is in 
perfect agreement with the one found before, i.e. the amplitude h~ is the same for the 
KA84 and VPI SP99 analysis. For di + (^2, ^3 and we encounter the same problems as 
discussed before and therefore refrain from giving the corresponding numbers. We note, 
however, that the values obtained for these LECs are very similar to the ones obtained 
for the Karlsruhe phase shifts, compare table |. For the reasons discussed in sec.|^, we 
consider the determination based on the KA84 phases more consistent (despite the fact 
that a part of the data set used in the Karlsruhe analysis is now outdated). As noted 
before, the dimension three LECs can mostly not be determined precisely. For achieving 
a better precision, one has to go to fourth order and fit in a larger region inside the 
Mandelstam triangle. 



5 Results and discussion 

In HBCHPT the quantities of interest are expressed in terms of the LECs. Explicit 
expressions for e.g. the subthreshold parameters and the cx-term (t(0) are given in §]. 
For convenience, we collect some of these expressions in the appendix. We concentrate here 
on a subset of quantities namely these which a) were found to be most problematic in the 
analysis using mostly data from the physical region and b) do only involve du — di^ from 
the dimension three LECs. P^ The data set eg . (p3D yields e.g. the following predictions 



for some of the low subthreshold coefficients (for the reasons mentioned above, we do not 
give any uncertainties): 

a+ = -1.32 M-^ , ato = 4.49 , a^^ = 0.97 , = 9.99 M'^ . (36) 

The subthreshold parameters are in good agreement with the Karlsruhe analysis, Qqq = 
-1.46 ± 0.10, afo = 4.66^3, a^i = 1.14 ± 0.02 and 600 = 10-36 ± 0.10, in appropriate 



units of the inverse pion mass. Of particular interest is the result for a^Q, which came 
out consistently too large in the chiral analysis based on the phase shifts, cf. table E.l in 
ref. Furthermore, the result for Bqq, which has not been given before in the literature, 
is in good agreement with the Karlsruhe analysis. While these results are very promising, 
one still has to check their stability by performing a fourth order calculation. 

Most striking, however, is the result for the pion-nucleon cr-term, 

a(0) = 40 MeV , (37) 

which agrees nicely with the dispersion theoretical analysis of ref.[^, a{0) = (44 ± 8 ± 
7) MeV, which is also based on the Karlsruhe phase shifts. On the other hand, the third 
order chiral perturbation theory analysis based on the phase shifts lead to much larger 
values of the sigma-term, a{0) = 59 MeV in ref. 0] and cr(0) = 70 MeV in ref.0 for the 

'^^'^Luckily, these two conditions are mutually consistent. 

'^^^Note that the Karlsruhe group does not give an uncertainty for this quantity. 
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corresponding central values of the LEG C\. This lends further credit to the statement that 
the chiral predictions for vrA^ scattering are most reliable inside the Mandelstam triangle 
since there the pertinent kinematical variables take their smallest possible values. For the 
reasons mentioned above, we find it difficult to determine a theoretical uncertainty to this 
number. We estimate the theoretical uncertainty (at this order) to be of the same size 
than in the analysis of ref.[^]. Using the formula derived in refs. p^, 25| 



(36±7)MeV , . ^iE^fW , (38) 

1 — y \'p\uu + dd\'p) 

we get for the strangeness content y = 0.10 ±0.17, somewhat smaller but compatible with 
the result of ref . |^ . If we insert the much larger value for ci obtained in the analysis based 
on the VPI phase shifts, we obtain naturally a much larger cx-term of about 209 MeV. 
This value is even larger than the one found in ref. , where the threshold nN phases from 
the VPI SP98 solution were analyzed. As we noted before, this particular LEG (together 
with C3) is sensitive to the subtraction function A~^{0,t). We also stated before that 
the subtraction function deduced from the SP99 solution leads to an unacceptably large 
deviation from the Adler consistency condition. Gonsequently, the very small isoscalar 



amplitude might not be well represented. Inserting such a large value for ci into eg. (|38|) , 
one obtains a huge strangeness content, y = 0.83 ± 0.04, which would mean that a large 
fraction of the nucleon mass is due to strange quark pairs. Such a picture of the nucleon 
is not tenable. However, one might argue that higher orders not taken into account 
in eq.(P5D might drastically alter the number (36 ± 7) MeV in the numerator of eq.(^) 
which is crucial in the link between a{0) and the strangeness content. While such a 
scenario is improbable, it can only be ruled out be a complete two-loop calculation of 
the scalar sector of baryon GHPT. Glearly, a full dispersive analysis using the modern 
data set is needed to further clarify this issue. We also mention that such a large cr-term 
would be at odds with all other information one has on nucleon matrix elements of the 
various strangeness operators like e.g. 57^ s or 57^75 s. The contribution of such matrix 
elements to observables like e.g. the proton magnetic moment seems to be fairly small. 



as indicated by recent measurements on parity-violation in electron scattering [|26|, 
What also should be done is to supply more stringent theoretical uncertainties. For 
that, one can not use the partial- wave analyses since these do not supply any error. 
Gonsequently, one has to reanalyze the pion-nucleon scattering data. This together with 
a more detailed discussion of the other amplitudes and subthreshold parameters will be 



given elsewhere [28 



6 Summary and conclusions 

In this paper, we have considered pion-nucleon scattering inside the Mandelstam plane 
from the point of dispersion and chiral perturbation theory. Dispersion relations are based 
on general principles like unitarity, crossing and analyticity and allow one to reconstruct 
the invariant amplitudes from the data, in particular also inside the Mandelstam triangle. 
We have performed such a calculation here based on the Karlsruhe partial-wave analysis 
as well as on the SP99 solution of the VPI group. On the other hand, chiral perturbation 
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theory is the effective field theory of the Standard Model at low energies and can be 
used to investigate the strictures from the spontaneous and explicit chiral symmetry 
violation. It is also based on general principles and admits a systematic power counting in 
small momenta and quark masses. Therefore, it should work best inside the Mandelstam 
plane. That this is indeed the case has been demonstrated here for heavy baryon chiral 
perturbation theory (after subtraction of the Born terms which need a special treatment 
because of the analytic structure). We have determined the so-called low-energy constants 
by a direct comparison of the amplitudes obtained from HBCHPT with the dispersive 
ones. From that, we can work out physical quantities and in particular, we find a pion- 
nucleon cx-term of a{0) = 40MeV, consistent with a previous dispersive analysis of the 
Karlsruhe data We have also found a much improved description of the subthreshold 
parameters a J) and fego- The VPI SP99 partial wave analysis leads to a much larger cx-term 
of about 200 MeV. This partial wave analysis is, however, less stringently constrained by 
strictures from analyticity, which we believe to be an essential ingredient for a precise 
determination of small quantities like the pion-nucleon cx-term. As noted before, the 
Karlsruhe analysis contains some data which are now considered inconsistent with the 
rest of the data set. Therefore, a full scale dispersive analysis based on the Karlsruhe 
method using only the accepted modern data set is called for. 
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A Expressions for some observables 

Here, we give the explicit expressions for the subthreshold parameters and the cx-term 
discussed in sect. |. These read 

< ^ ?^(e,-2.).||!(.^,.?).0,M^). ,A.l) 

cx(O) = -^c^M'--^ + 0{M'). (A.5) 
Note that we have corrected for a typographical error that appeared in eq.(E.5) of ref.[S]. 
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